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Microfluidic devices can be used to produce highly controlled and monodisperse dou- 
ble or multiple emulsions. The presence of inner drops inside a jet of the middle phase 
introduces deformations in the jet, which leads to breakup into monodisperse double 
emulsions. However, the ability to generate double emulsions can be compromised 
when the interfacial tension between the middle and outer phases is low, leading to 
flow with high capillary and Weber numbers. In this case, the interface between the 
fluids is initially deformed by the inner drops but the jet does not break into drops. 
Instead, the jet becomes highly corrugated, which prevents formation of controlled 
double emulsions. We show using numerical calculations that the corrugations are 
caused by the inner drops perturbing the interface and the perturbations are then 
advected by the flow into complex shapes. 



^-'Present adress: School of Physics, Georgia Institute of Technology, Atlanta, Georgia 30332, USA. 
Present adress: Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, 

New Jersey 08544, USA. 
Electronic mail: 



weitz@seas.harvard.edu 



1 



I. INTRODUCTION 



Multiphase flows are widely used industrially for a plethora of applications including oil 
extraction, formulation of food, and home and personal care products. A good understanding 
of multiphase flows is needed to facilitate the use of these highly complex systems in the 
various applications. For example, the use of core-annular flows facilitates water-lubricated 
transport of heavy viscous oils by providing new approaches to mitigate problems such as 
fouling of oil pipes.'^ Depending on the properties of the flow, such as pressure gradients 
and fluid velocities, a variety of flow structures including corrugated interfaces has been 
observed. While the study of jets and drops in multiphase systems has facilitated numerous 
conventional applications,'^ the emergence of microfluidics creates new opportunities to study 
multiphase flow in microscale systems.!^ 

Microfluidic methods for manipulating multiphase, water-and-oil systems have a great 
potential for applications involving encapsulation. One particularly promising way to ac- 
complish such encapsulation is through the controlled generation of multiple emulsions,'^^ 
which can be done in various conflned microfluidic channels with uniform depths, such as 
T-junction^ and flow-focusing channels.'^'^ Another approach to form these emulsion is to 
direct droplet breakup in capillaries aligned axisymmetrically.'^ The breakup of the inner 
phase and the middle phase can either happen simultaneously or in two steps. In the latter 
case, the inner drops are flrst generated and trigger the breakup of the middle jet into double 
emulsions due to interfacial tension effects. Tapered capillaries are typically used in these 
devices so that droplets slow down after generation, resulting in lower viscous stresses, which 
help prevent destabilization of the droplets due to these stresses. The microfluidic technique 
allows the formation of monodisperse double emulsion with a high degree of control over 
the size and shape of the emulsion drops, and with excellent encapsulation efficiency.'^^'^ 
Thus it provides a versatile platform for studying the dynamics of drop formatiod ^^ * ^^ ! and 
for fabricating novel functional materials .'^^'^^'^ 

Despite the high degree of control afforded by the technique, the multiple emulsions 
formed are typically not stable due to coalescence between miscible phases. For example, in 
a double emulsion, the innermost drop can coalesce with the outermost, continuous-phase 
fluid, destroying the multiple emulsion. A common strategy to improve emulsion stability 
is to introduce surfactants, which generally help to suppress coalescence. However, the in- 



troduction of surfactants also leads to lower interfacial tensions; the dynamics associated 
with these lower interfacial tensions can result in a variety of fluid flow patterns and in- 
stabihties that disrupt the regular formation of the double emulsion. Consequently, these 
effects can severely affect the encapsulation efficiency and prevent prolonged operation of 
the processes, which makes microfiuidic emulsification unfeasible for industrial applications. 
It is, therefore, important to understand the origin and mechanisms of formation of these 
flow patterns. 

In this paper, we study the corrugations observed along the fluid-fluid interface inside a 
glass capillary device when the interfacial tension between the middle phase and the outer 
phase is very low, due to the introduction of surfactants. We operate the device in a regime 
where the innermost fluid forms drops by dripping right at the oriflce, whereas the middle 
fluid forms a jet. We study experimentally the interfacial corrugations, which are initially 
triggered by advection of the innermost drops inside the jet. Due to the low interfacial 
tension between the middle and outer phases, the perturbations of the interface can be 
advected into complex shapes by the flow. We observe the phenomenon when the surface 
tension is low compared to both inertial and viscous forces; as a result, both the capillary and 
Weber numbers are very large. We suggest an explanation of the mechanism of formation 
and evolution of the corrugations, which is validated by reproducing the corrugations in a 
numerical model that agrees qualitatively and quantitatively with the experiments. 



II. EXPERIMENTAL 

All experiments were carried out in glass microcapillary devices,'^ shown schematically 
in Fig. [l](a). The circular capillaries, with inner and outer diameters of 0.58 and 1.0 mm 
(World Precision Instruments, Inc.), were tapered to the desired shapes with a micropipette 
puller (P-97, Sutter Instrument, Inc.) and a microforge (Narishige International USA, Inc.). 
Two tapered capillaries with different tip diameters were aligned inside a square glass cap- 
illary (Altantic International Technology, Inc.) with an inner dimension of 1.05 mm. The 
innermost dispersed phase flows through the round capillary with the smaller-diameter tip; 
this is the injection tube. The second round capillary, with the larger-diameter tip collects 
the double emulsion and directs it to the exit of the device. The diameter of the collection 
tube ranged from 80 to 200 /xm, while the diameter Di of the injection tube ranged from 40 
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to 100 ^m. A transparent epoxy resin (5 Minute Epoxy, Devcon) was used to seal the capil- 
laries where necessary. Positive syringe pumps (PHD 2000 series, Harvard Apparatus) were 
used to deliver the different phases at the desired ffow rates. The ffow behavior inside the 
devices was monitored with a lOx objective using an inverted microscope (Leica, DMIRBE) 
equipped with a high-speed camera (Phantom, V5, V7, or V9). Unless otherwise noted, the 
inner, middle, and outer phases consisted of water, 2% — 10% (w/w) Span 80 surfactant in 
dodecane, and water with — 8 mM sodium dodecyl sulfate (SDS), respectively. Solutions 
and solvents were all filtered before introduction into the glass microcapillary devices. 

Unless otherwise stated, all chemicals were purchased from Sigma-Aldrich. Water with 
a resistivity of 18.2 MQ cm~^ was acquired from a Millipore Milli-Q system. Whenever 
possible, characteristic interfacial tensions were measured by forming a pendant drop of the 
denser phase at the tip of a blunt stainless steel needle (McMaster-Carr, 20 Gauge) immersed 
in the second phase and performing drop-shape analysis.!^ Low interfacial tensions that could 
not be precisely determined using the pendant drop method were measured using a spinning 
drop tensiometer (SITE 100, Kruss GmbH). Viscosities were measured with a Ubbelohde 
capillary viscometer mounted in a temperature-controlled water bath. 

III. RESULTS AND DISCUSSION 

In a typical system for generation of emulsions, the inner phase is typically at a low flow 
rate, allowing droplets of the inner phase to be formed at the orifice of the injection tube in 
a flow regime known as dripping. By contrast, the middle phase is typically injected at high 
flow rates; as a result, a jet of the middle phase is formed inside the collection tube; thus 
it is operated in a regime known as jetting. The advection of the inner droplets inside the 
middle jet causes the interface between the middle and the outer phases to deform. With a 
typical interfacial tension, for example, 7 = 2.7 ± 0.2 mN m~^, for the experiment shown in 
Fig. [l](b), the deformation grows and eventually causes the middle jet to break up, thereby 
forming double emulsion drops. The breakup of the middle jet occurs as a consequence of 
the interfacial tension and viscous stresses due to the flow of the continuous phase. While 
interfacial tension tends to favor droplet formation, viscous stresses induce the interface to 
follow the flow proflle inside the channel. The balance of these two effects is characterized 
by the capillary number of the continuous phase: Ca = rjoutVout/l, where rjout and Vout 
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FIG. 1. Double emulsions and the formation of corrugated interfaces, (a) Schematic of the exper- 
imental setup, (b) Generation of monodisperse double emulsion droplets. The inner, middle, and 
outer fluids are, respectively, water, 10% (w/w) Span 80 in dodecane, and water at flow rates of 
1000, 1500, and 5000 h~^, respectively. In this system, the interfacial tension between the mid- 
dle and outer phases is 7 = 2.7 it 0.2 mN m~^. Re ~ 15, Ca ~ 0.15, and We ~ 0.6. (c) Formation 
of corrugations when the interfacial tension between the middle and outer phases is significantly 
reduced. The inner, middle, and outer fluids are water, 2% (w/w) Span 80 in dodecane and water 
with 4 mM SDS and the flow rates are 500, 10^, and 10^ h^^ for the inner, middle, and outer 
fluids, respectively. In this system, the interfacial tension between the middle and outer phases is 
30 lb 4 fiN m^^, Re ~ 38, Ca ~ 37, and We ~ 600. The radius of the capillary at its entrance Rq 
is 38.8 ^m. Scale bars are 100 //m. 

are the viscosity and velocity of the outer continuous phase, respectively. The flow in the 
experiment shown in Fig. has a low capillary number, Ca ~ 0.15. Apart from viscous 
stresses, inertial effects can also compete with surface tension to prevent droplet breakup. If 
the jet moves faster than the continuous phase, the velocity difference induces a large shear 
at the interface that decelerates the jet, causing it to widen. The balance of this inertial 
effect and surface tension is characterized by the Weber number of the jet-forming phase, 
which includes both the inner and middle phases: We = pm dVin^ /'j, where pm and Vm are 
the density and velocity of the jet-forming phase, respectively. The flow in the experiment 
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shown in Fig. [T](b) has a low Weber number, We ~ 0.6. 

However, the performance of the glass microcapillary devices in generating controlled 
double emulsions is compromised when the interfacial tension between the middle and the 
outer phases is low such that Ca, We > 0(1). Under such flow conditions, the middle jet 
does not break up; instead a highly corrugated interface is formed, as shown in Fig. ^c). 
Here the interfacial tension is 30 ± 4 /iN m~^, Ca ~ 37, and We ~ 600; both Ca and We 
are about 100 times the values for the flow shown in Fig. [l](b). 

To understand the formation of the highly corrugated interface, we examine the onset 
of the corrugation with optical microscopy. Since the interfacial tension between the inner- 
most droplet phase and the middle phase is about two orders of magnitude higher than that 
between the middle phase and the outermost continuous phase, the interface between the 
middle and outermost phases is significantly more deformable than that between the inner- 
most and middle phases. As a result, when an inner droplet is about to pinch off from the 
injection tube in the dripping regime, the interface between the middle and the outer phases 
starts to become deformed. At this point, the inner droplet is still attached to the tip of the 
injection tube; however, the droplet size is already larger than the size of an unperturbed 
middle jet. This geometry causes the middle jet to curve around the inner droplet, as shown 
in Fig. ga). 

After the complete pinch-off, the inner droplet continues to flow downstream at a higher 
velocity than the deformed interface, which lags behind, as shown in Fig. ^h). As the fluids 
flow downstream, the shape of the deformed interface evolves such that the front of the 
deformation has a smaller diameter than its back [Fig. [2|c)]. We also measure the number 
of corrugations as a function of the number of inner drops generated, as shown in Fig. [3] The 
frequency of the corrugations approximately matches the generation frequency of the inner 
droplets. These observations suggest that the corrugations are initiated by the formation of 
the inner drops, rather than simply arising as an intrinsic fluid or interfacial instability such 
as the Rayleigh capillary instability. While the formation of the corrugations is triggered by 
the motion of the inner drops, the high curvature in the resultant highly corrugated jet is 
sustained because of the ultralow interfacial tension of the interface. 

Two potential driving forces for the evolution of the interface to its final shape are the 
disturbance flow due to the moving inner drop and advection of the initially deformed 
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FIG. 2. (a) Formation of an inner droplet in the dripping regime. The dotted white hne indicates 
the location of the injection tube. The scale bar is 100 /im and it applies to [(a)- (c)]. (b) Deforma- 
tion of the middle jet by the pinch-off of the inner droplet, (c) Evolution of the deformed middle jet 
into a corrugation, (d) Proposed mechanism of the corrugation formation. The interface between 
the middle and outer phases is initially deformed by the droplet and then evolves following the 
parabolic flow configuration due to the flow in a widening channel. 

interface in a pressure-driven Poiseuille flow in a widening capillary. If the former effect 
dominates, the inner droplet should continue to significantly deform the interface between the 
middle and outer fluids after the drop formation. However, the magnitude of the deformation 
caused directly by the drop is inconsequential when compared to that of the corrugation, 
as suggested by the largely undeformed interface near the inner droplet in Fig. ^c). Thus, 
it seems unlikely that the disturbance flow due to the moving inner droplet is the main 
driving force. Therefore, we hypothesize that the final shape of the corrugations is caused 
by advection of an initially perturbed interface in a Poiseuille flow in a widening channel, 
as illustrated schematically in Fig. ^d). 

Based on these observations, we suggest a mechanism for the formation and evolution of 
the corrugations: the inner fluid breaks up into droplets and due to the incompressibility of 
the middle fluid, the interface between the middle and outer fluids is deformed into a nearly 
spherical shape around the sphere. Because of the low interfacial tension, the viscous stresses 
and inertial forces dominate. Thus the deformed interface passively flows downstream ac- 
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FIG. 3. Number of corrugations observed, N corrugation-, as a function of the number of inner drops, 
Niirop, for a given length of the device. The sohd hne has a slope of unity. Since the inner drops flow 
at a higher speed than the corrugations and the overall speed decreases as they flow downstream in 
the widening channel, the number of corrugations counted in a typical image is a little larger than 
that of the inner drops. Each data point was obtained by counting the number of corrugations and 
the number of inner drops in at least five images. Inset: Snapshots of the flow when the flow rates 
of the inner phase are (a) 0, (b) 5, (c) 50, and (d) 500 fiL h~^ while keeping the flow rates of the 
middle and outer phases at 10^ /nL h~^ and 1.5 x 10^ /uL h~^, respectively. Scale bars are 100 /xm. 

cording to the two-phase Poiseuille flow proflle. Then, due to the approximately parabolic 
velocity distribution, the part of the curved interface that is closer to the channel wall flows 
at a lower velocity than the part near axis of the channel. As a result, the deformed interface 
develops into arrowlike, folded corrugations, as shown in Fig. ^c). 

To conflrm this proposed mechanism, we construct a model for the occurrence of the 
interfacial corrugations of the experimental system that highlights the most important ideas. 
In our experimental system, the viscosities of the middle and the outer phases are 1.44 and 
1.0 cP, respectively. Since the viscosity ratio is close to one, we assume a single-phase 
pressure-driven velocity fleld in our numerical calculations; this simplifles the calculations 
signiflcantly while still capturing the main features of the experimental observations. The 
initial condition consists of two lines of imaginary tracer particles that approximate the 
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initial shape of the experimentally observed interface immediately after pinch-off of the 
inner droplet, as shown in Figs. [2]^a) and|4|^a). The only role that the inner drop plays in 
our calculations is to initially perturb the interface. The tracers are then advected by the 
parabolic velocity distribution in the widening channel 



where R{z) = az + Ro is the inner radius of the capillary at the axial distance z from the 
entrance of the capillary and r is the radial distance of the tracer from the centerline. The 
capillary has a radius of Rq at its entrance, and widens with a slope of a, which has a typical 
value of about 0.23, as measured experimentally. 

As the perturbed interface flows downstream, the tracers that are furthest from the 
capillary axis lag behind those that are closer to the capillary axis. As a result, the initially 
spherical deformation sharpens, as shown in Fig. |4](b). The part of the interface further 
away from the axis continues to lag behind and the corrugation folds, as shown in Fig. ^c). 
The three-dimensional representation of the shape based on the results from our numerical 
calculations closely resembles the experimentally observed shape, as shown in Figs, ^d) 



We also track the positions of the tip of a corrugation obtained from the numerical 
calculations over time, neglecting the effects of interfacial tension. For a comparison with 
the experimental results, we normalize both the experimentally obtained and numerically 
calculated data. The x-positions of the tip are normalized by the radius of the capillary at 
the entrance, Rq. The timescales are normalized by a reference timescale, {Ro/ Umax, sphere), 
where Umax,sphere IS the speed of the sphere at the entrance of the channel. We then plot the 
normalized positions of an inner drop and the tip of a corrugation as a function of normalized 
time in Fig. |5} The experimental observed positions of the tip, represented by the square 
symbols, agree very well with the numerically calculated positions of the tip, shown with a 
solid line in the plot. This result supports our hypothesis that the mechanism that produces 
the corrugations is folding of an initially perturbed shape by the slowly varying parabolic 
velocity profile. The result of the calculations deviates slightly from experiments at long 
time, as shown by the increasing gap between the symbols and the solid line in Fig. |5| 




(1) 



andge). 
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FIG. 4. (a) The initial shape of the interface after pinch-off of an inner droplet. The shape is used 
as the initial condition for the numerical calculation, (b) Shape of the interface developed under 
pressure-driven flow in the widening channel, (c) Shape of the interface further downstream from 
(b). (d) Three-dimensional shape of the interface by rotating the two-dimensional shape in (c) by 
360° . (e) Image of the corrugations observed experimentally. The shape of the corrugations is 
captured qualitatively by the calculation, as shown in (d). Scale bar is 100 //m. 

The deviation likely results from effects on the flow due of the motion of the inner droplets, 
which have not been included in our simple model. 

To test the robustness of the proposed mechanism, we also numerically calculate the 
shape of the perturbed interface obtained in a device where the injection and collection 
tubes are slightly off-center. This feature is often observed experimentally and leads to 
axial asymmetry in the droplet generation process. As a result of the misalignment, the 
inner drops are closer to one side of the collection tube and the initial shape of the interface 
is changed accordingly. To confirm the validity of our model, we perform our numerical 
calculations with an initial condition where the two lines of imaginary tracer particles are 
asymmetrically positioned to approximate the initial shape of the interface deformed by an 
off-center inner drop. From the results of our calculation, we obtain corrugations that are 
axially asymmetric about the capillary axis, as shown in Fig. [6|^a). The shape obtained 
is similar to the experimentally observed asymmetric corrugations shown in Fig. [6]^b). To 
highlight the close resemblance, the boxes of the same size and shape, drawn with black dot- 
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FIG. 5. Measured positions of the inner drop (Q) and the tip of the corrugation (□) as a function 
of normahzed time, t/ [Rq/ Umax, sphere)- The sohd Hne is the prediction for the position of the tip 
of the corrugation from the numerical calculation with no fitting parameters. Inset: a snapshot of 
the corrugated interface, with the positions of the inner drop (Q) and the tip of the corrugation 
(□) marked. The radius of the capillary widens from Rq at the entrance at a slope of about 0.23. 
Scale bar is 100 /im. 

ted lines, are placed around both the experimentally observed and the numerically calculated 
corrugations in Figs. [6|^a) and[6|^b). 

To summarize, we have observed corrugated interfaces in a multiphase core- annular flow 
at low interfacial tension in glass capillary microfluidic devices designed for the generation 
of double emulsions. We present qualitative evidence that the initial perturbations of the 
compound jet are introduced by the advection of the small inner drops inside the jet. As a 
result of the low interfacial tension between the middle and outer phases, the perturbations 
are advected into complex shapes by the flow; the high curvature of the resultant highly cor- 
rugated shape is not suppressed since the Weber numbers are large (low values of the surface 
tension). We also perform numerical calculations in which a single-phase pressure-driven ve- 
locity field advects tracers that outline the interface. Despite the simplifying assumption, 
our calculations capture the qualitative and quantitative behaviors of the experimentally 
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FIG. 6. (a) Numerical calculations of the corrugations on the surface of a misaligned middle 
jet. (b) Experimentally observed corrugations on the surface of a slightly off-centered middle jet. 
Misalignments of the capillaries in the device lead to axial asymmetry in the positions of the 
different phases, often resulting in asymmetric corrugations. Such asymmetry is captured by the 
calculation. To highlight the close resemblance, a box of the same size and shape drawn with black 
dotted lines are placed around both the experimentally observed and the numerically calculated 
corrugations in A and B. Scale bar is 100 /xm. 

observed corrugations. The improved understanding of this phenomenon will help provide 
a higher degree of control over multiple emulsions generated with microfluidic techniques. 
Moreover, our study also suggests a new route to create highly corrugated interfaces. With 
an appropriate strategy for solidifying the middle jet, our system may lead to fibers with 
significantly enhanced surface area for biomedical and catalytic applications that requires 
large surfaces for improved activities. 
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